Chaos-assisted, broadband trapping of light in optical resonators 
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Chaos is a phenomenon that occurs in many aspects of contemporary science. In classical dy- 
namics, chaos is defined as a hypersensitivity to initial conditions. The presence of chaos is often 
unwanted, as it introduces unpredictability, which makes it difficult to predict or explain exper- 
imental results. Conversely, we demonstrate here how chaos can be used to enhance the ability 
of an optical resonator to store energy. We combine analytic theory with ab-initio simulations and 
experiments in photonic crystal resonators to show that a chaotic resonator can store six times more 
energy than its classical counterpart of the same volume. We explain the observed increase with 
the equipartition of energy among all degrees of freedom of the chaotic resonator, i.e. the cavity 
modes, which is evident from the convergence of their lifetime towards a single value. A compelling 
illustration of the theory is provided by demonstrating enhanced absorption in deformed polystyrene 
microspheres. 



Introduction 

The enhanced interaction between light and matter in 
optical cavity resonators is a subject of a great interest 
that impacts on many areas of condensed matter physics 
such as cavity electrodynamics [1], quantum and non- 
linear optics, but also on more applied aspects such as 
optical signal processing 0-11] and resonantly enhanced 
optical absorption [B'l . All these applications are enabled 
by highly optimized optical resonators that can efficiently 
trap electromagnetic energy in narrow frequency bands. 
In conventional cavities, there is a simple trade-off be- 
tween bandwidth and the enhancement of trapped en- 
ergy; the higher the enhancement, the narrower the band- 
width. A great challenge in the field is therefore to de- 
velop a new generation of cavities that are able to break 
this fixed relationship and to store more energy in a given 
bandwidth window than conventional cavities would al- 
low. 

The maximum power that can be transferred into a con- 
ventional resonator depends on the coupling coefficient 
and loss of each given mode and tends to vary across 
the mode spectrum, especially when broadband opera- 
tion (A A « lOO's of nm) is considered. In addition, clas- 
sical 2D and 3D geometries tend to accommodate modes 
with very different lifetimes in the same spectral region, 
a good example being the widely used photonic crystal 
"L3" type cavity that features modes of very different 
Q-factor closely spaced in frequency @. Therefore, the 
use of classical resonators for broadband energy storage 
is limited. Here, we overcome this intrinsic limitation 
by exploiting specific shape deformations that support 
chaotic trajectories for light rays. Surprisingly, we note 
that in a chaotic cavity, the lifetimes of all the modes tend 
towards a common value, thus improving the transfer of 
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energy into the cavity and increasing the energy-storage 
capability of the cavity. Such chaotic resonators 01 have 
been well exploited in the field of laser devices |8l-[l5|. 
We also note that the mode spectrum of deformed mi- 
crosphere resonators has been studied recently with re- 
markable changes in the Q-factors being observed 

PI- 

However, despite this large body of literature, nothing is 
known about the capacity of such resonators to store and 
collect light energy over a broad spectral range. 
The increased energy storage capacity of a chaotic res- 
onator compared to a classical one can be intuitively ex- 
plained by adopting a ray optics approach and consider- 
ing that a suitable shape deformation is accompanied by 
the breaking of symmetry in the structure. As a con- 
sequence, the deformed resonator cannot support any 
cyclic motion of light, thus the trajectory of light rays 
changes from regular to random, statistically resulting 
in a larger lifetime of the photons in the cavity [7]. In 
order to clarify this result, we start our analysis from a 
symmetric (classical) resonator and observe its capacity 
to trap energy as it is deformed. We define a deforma- 
tion parameter, a, which we will use in the following 
description as a handle to deform any given geometry, 
with a — indicating the original, undeformed structure 
and a > indicating a proportionally deformed geome- 
try. For example, a circle would be described by a = 0, 
and a deformed circle by a > 0, with a describing the 
degree of deformation. This parameter is generic and can 
be used to describe the deformation of any type of res- 
onator, e.g. square or disk in 2D or cube or sphere in 3D. 
In general, the larger a, the larger the degree of chaos, 
until a saturation value is reached. The parameter a is 
defined in Eq. ([1]) below. 

We begin our analysis with a single ab-initio numerical 
experiment, which demonstrates a six-fold increase of the 
energy stored inside a chaotic resonator with respect to 
a classical counterpart of the same volume. Such an an- 
alytic treatment allows us to study the lifetimes of the 
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FIG. 1. (a) Snapshot of the electromagnetic energy density 
H distribution after t — 45/s in a resonator (dashed hne) with 
V — SO/im^ for (a) a = and (b) a = 1.5; (c) time evolution 
of the average electromagnetic energy {£) for different a. 



electromagnetic modes excited in the cavity and to exam- 
ine the coherent buildup of energy inside the system. We 
find that energy gets uniformly distributed among the 
spectral degrees of freedom and that the cavity is able to 
store the same amount of energy for a given wavelength 
interval. We support our theoretical results with a set 
of experiments both in 2D and 3D geometries. Two- 
dimensional resonators were fabricated in planar pho- 
tonic crystals (PhCs) and analyzed by pump and probe 
transmission measurements. We observed and character- 
ized the resulting modal lifetimes for varying degrees of 
chaos. While the PhC environment provides fine con- 
trol over the shape of the resonator, hence its chaotic 
behavior, it is difficult to measure the energy stored in- 
side the cavity directly. For this reason, we also studied 
a three-dimensional arrangement of polystyrene spheres, 
each suitably modified to exhibit chaotic behavior. De- 
spite the high transparency of polystyrene, absorption 
measurements showed a broadband absorption enhance- 
ment of ~ 10% across its entire absorption bandwidth 
(« 450nm), due to the larger energy trapped in the 
chaotically deformed resonator. It is worthwhile to stress 
that the only reason for the increased amount of energy 
trapped in the resonator is the onset of chaos caused by 
the deformation of the spheres, with the amount of ma- 
terial being kept strictly constant. 



A single ab-initio experiment: broadband light 
confinement and energy equipartition 

The theory is more easily developed in 2D starting 
from a circular resonator. Our setup consists of a silicon 
dielectric resonator with air holes , whose shape is defined 
by the following function in polar coordinates {pc,Oc)- 



Pc=\ -+acos6lc, 0<6l, <27r, (1) 

V TT Z 



with A the resonator area and a > the single parameter 
that controls the resonator shape. Equation ([T]) belongs 
to the family of analytic curves investigated by Robnik: 
for a > 0, Eq. ([T|) supports chaos in the trajectory of 
light rays, which randomly bounce inside the resonator 
pTj . From a physical perspective, the shape defined by 
Eq. IT]) is equivalent to an asymmetric deformation of 
a disk, and can be realized experimentally with conven- 
tional nanofabrication tools. 

In our simulations, we fixed the resonator area to A = 
30 firri^ and numerically calculated the electromagnetic 
energy £ stored inside the resonator for varying values of 
a. Although the value of a in Eq. ([T]) is not bound to an 
upper value, we note that the system reaches saturation 
for a maximum degree of chaos described by a = 1.5 (see 
Materials). We therefore restrict our numerical analysis 
to a e [0,1.5]. The calculation of the electromagnetic 
energy £{t) — Jj, Hdp is performed by a numerical inte- 
gration of the energy density H(/9; t) = i(E-D + H-B) in 
the volume V encompassed by the resonator and defined 
by p < Pc and < < 27r. We simulated the input from a 
supercontinuum source in the wavelength range between 
A — SOOnrn and A = ISOOnm, which simulates a broad- 
band source such as sunlight. Figure [T^-b shows a time 
snapshot of the spatial distribution of after t = 45 fs, 
and illustrates how the geometry of the resonator changes 
with a (Fig. [T^-b, dashed line). As seen, a small defor- 
mation in the shape (at a constant volume) already yields 
a radically different behavior in the distribution of light 
energy inside the resonator. Figure [T];, conversely, dis- 
plays the time-averaged energy {£) — j dt'£{t') evolu- 
tion with increasing parameter a. Quite remarkably, the 
introduction of chaos into the motion of light is accom- 
panied by a dramatic change of the energy stored inside 
the resonator, whose steady state regime — when the in- 
sertion of energy balances radiation losses — shows an 
6-fold increase already for a — 1.5 (Fig. ^p). This energy 
accumulation grows with the deformation, and monoton- 
ically rises as a is increased from a = 0toQ! = 1.5. 

We use Time Dependent Coupled Mode Theory (TD- 
CMT) [l^ to develop a simple and reliable model for the 
light-resonator interaction. The system can be modeled 
as a side coupled resonator, whose dynamical equations 
can be easily solved to get the energy £k and the power 
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FIG. 2. (a) Log-plot of the energy £ relaxation dynamics for a = and a = 1.25: FDTD results (markers) and Prony 
exponential fits (solid lines) ; (b) normalized distribution of the difference between the maximum and minimum decay constants 
A(q:)/A(0) = A/Ao versus a; (c) steady state energy distribution {£s} = {£){t = 45/s) for varying values of a. The insets 
in (c) sketch the effects of the convergence of the decay constants Tk on the power Vk transferred into the structure, whose 
behavior versus Te/rfco is plotted as a solid line; in the non-chaotic situation (left) the distribution of Te/rko (left, colored 
arrows) is very broad and only a few frequencies are able to efficiently transfer power; in the strongly chaotic case, conversely, 
a condensed distribution of Te/rko yields the same contribution for all wavelengths and a much larger number of frequencies 
can contribute to storing energy in the structure. 



Vk stored into the the k—th mode (see Materials): 

= 'P^= .^%2 \S\', (2) 

with \S\ the input source power, = ^ -I- the mode 
decay rate, with l/rko the intrinsic cavity decay rate of 
the k-th mode and I/te the escape rate due to coupling 
with the environment. The power Vk strongly depends 
on the ratio between the radiation and the coupling loss 
through the parameter r^o / T'e , achieving the maximum 
value of Vk = 0.5|5p when Tg = Tko- Outside this match- 
ing condition, the power coupled into the structure de- 
creases very fast. To study how the dynamics of the decay 
rates can be affected by chaos, we began by calculating 
the evolution of the decay rates Tko for different a in the 
resonator defined by Eq. ([T]) . In a series of FDTD simu- 
lations, we first excited the resonator with a source and 
then monitored the energy evolution £{t) when the source 
was switched off. The decay constants Tko are extracted 
from the time energy evolution SU) = l^fcj^e"^*/'^'''' 
by applying the Prony method Fig. [5^ shows two 
typical examples of numerical simulations (markers) and 
fit (continuous line) for two different values of a. In gen- 
eral, modes of different frequency exhibit different de- 
cay rates I/t/cq. However, when strong chaos is gener- 
ated in the structure, the distribution of t^o converges 
towards a frequency independent delta function. In fact, 
calculations which start from initial conditions that al- 
ready belong to the chaotic region of phase space will 
always be distributed according to the same probability 
distribution, which defines the so called natural invari- 
ant measures of the system [20| . As a result, when the 
entire phase-space is dominated by chaos, we observe the 
same evolution of the decay rates for all possible initial 
conditions towards a frequency-independent decay rate 
Tko — tq. a convenient way of highlighting this dynam- 
ics is to plot the difference between the maximum and the 



minimum decay constant: A(a) = max(Tfeo) — niin(Tfeo), 
for different values of a (Fig. [^b)- We clearly observe a 
transition scenario: below a*, the dynamics shows sim- 
ply an oscillation of A(a) around the same average value, 
while above the chaos threshold a > a* , a, clear conver- 
gence of A — is observed. The value a* depends on 
the specific geometry of the chaotic resonator, and can 
be assessed by calculating the relative area of the system 
phase space that encompasses chaos (see Materials) . The 
effect of this convergence towards a single lifetime of all 
the modes on the energy collected by the resonator can 
be readily evaluated from Eq. ^ . The total power trans- 
ferred into the structure, in particular, then becomes fre- 
quency independent Vk ^ Vo = 2|S'P(l + Te/to) Te/ro, 
and every mode contributes to the same extent to stor- 
ing energy inside the resonator. When many modes are 
present in the resonator, their large number results in a 
coherent buildup process that leads to a significant ac- 
cumulation of energy (Fig. [T}d and Fig. [2}:). In the 
non-chaotic case, conversely, much fewer modes are able 
to efficiently transfer energy into the resonator due to 
the mismatch between Tko and Tg (Fig. inset) and 
the system can store relatively less energy (Fig. [TJd and 
Fig. Hi). 

The chaos-assisted energy buildup process observed 
when a > a* originates from the fundamental thermo- 
dynamic principle of equipartition, which can be high- 
lighted using Eq. By substituting r^o — tq into the 
left hand of equations and assuming a dense distri- 
bution of modes, with wavelength separation Afe+i ~Xk — 
dX <C A, we obtain: 

^ = const. = So, (3) 

which can be regarded as an equipartition theorem, 
with the energy equally distributed among all degrees 
of freedom — i.e., the spectral wavelengths — due to the 
strongly chaotic nature of the system. Equipartition is 



4 







6 












(b) 


A?.=875nm 










AI=625nm 






/ AI=250nm 








c 










) 


10 20 , 
time [ps] 


y ■ 










theory 










■ FDTD, a=l,5 










• FDTD, a=0 




• 


• • 


• 





AX [|am] 



FIG. 3. Results for a variable bandwidth source: (a) FDTD 
calculated average energy (£") versus time for a = 1.5 and 
different normalized bandwidth AA symmetrically centered at 
800nm; (b) FDTD computed steady state average energy (fs) 
versus bandwidth AA for a = (circle markers) and a = 1.5 
(square markers). In (a) the solid line indicates the behavior 
predicted by Eq. The energy (fs) is normalized to the 

maximum value Em. attained for ol = 1.5 and AA = 1/xm. 



at the foundation of classical statistical mechanics, and 
forms the basis for thermodynamic ensembles and the 
observation of different phases of matter [2l[. Applied 
to Photonics, we have the remarkable opportunity of ex- 
ploiting this principle for enhancing the energy confine- 
ment properties of photonic structures. As already dis- 
cussed above, this opportunity was confirmed by calcu- 
lating the energy stored for a broadband source in the 
two limiting conditions of a fully chaotic geometry (Fig. 
[5^ square markers) and a non-chaotic geometry (Fig. [3^ 
circle markers) , and finding a six- fold enhancement in the 
chaotic case. 



2D experiments: Controlling the convergence 
towards a single modal lifetime 

Energy equipartition and chaotic energy harvesting 
both originate in the convergence of the modal lifetimes 
towards a single value r^o . We therefore begin our exper- 
imental campaign by investigating the occurrence of this 
phenomenon in real structures. We designed and exper- 
imentally realized a series of planar 2D stadium-shaped 
resonators in planar photonic crystals (PhCs) fabricated 
in Silicon on Insulator (SOI). The substrate consists of a 
220nm thick silicon layer on a 2/im thick insulator buried 
oxide. The patterns were written into ZEP resist on a 
modified LEO/RAITH system with 2 nm step-size and 
etched with a 50:50 mixture of SF6 and CIIF3 gases in a 
reactive ion etching machine. After stripping the residual 
resist, the sample was cleaved for end-fire coupling. Fig- 
ures (|l])a-c show a set of SEM micrographs and illustrate 



how the shape of the resonator evolves with a change 
of the deformation parameter a. The shape starts as a 
regular square (Fig. for a=0 and turns into a fully 
chaotic stadium-shaped resonator (Fig. HJ:) for a=1.5, 
with strong chaos alrea dy developed for a > a* = 0.5. 
The parameter a — y/ A/{A + no?) guarantees a con- 
stant resonator area A as a is varied (Fig. HJa). As the 
PhC lattice was designed to exhibit a bandgap around 
1.5^m with a bandwidth of « 400nm, any electromag- 
netic wave in this range is perfectly reflected at the PhC 
boundaries, and can only escape via the input /output 
waveguides of the structure or by scattering at imperfec- 
tions. The system can therefore be described as a 2D 
resonator with measurable losses and the modal decay 
constants can be extracted from the transmitted spec- 
trum. We employed the multiscale analysis described 
in [22], which provides an excellent reconstruction tech- 
nique even when the resonances overlap. The method 
fits the spectrum via a sum of suitable wavelet func- 
tions, allowing to compute the wavelength position Aq 
and the width i5A of each mode. The resonance widths 
8\ are then inversely proportional to the modal decay 
rates l/rfe = (5w = cSX/tlXq [l^. In order to collect 
statistically relevant data, we fabricated several samples 
with different areas: A = 400^^2, A = SOO/im^ and 
A ~ 1200/1771^, and for each area we considered five dif- 
ferent degrees of chaos, expressed via the deformation 
parameters a — a — 0.25, a — 0.5, a — 0.75 and 
a = 1. To characterize the samples, we used polarized 
light in the wavelength range 1520 nm to 1620 nm. Fig- 
ure [JJi shows a portion of one of the measured spectra, 
together with its reconstruction via multiscale analysis 
[231, which shows a perfect reproduction of the experi- 
mental results. All spectra (not shown here) have been 
reconstructed with the same level of accuracy. We are 
able to extract « 2000 resonances for each a , which al- 
lows us to extract statistically relevant trends. Figure HJj 
displays the resonance linewidth probability distribution 
P{5uj) for a = 1 (stadium) and a = (square). We note 
that for a = 1, the resonances are strongly converging 
towards a value of « 5 • 10~'^ THz with only negligi- 
ble contributions arising from the short lived modes that 
are characterized by a larger 5X. Conversely, in the non 
chaotic regime, we observed the presence of many short 
lived resonances, indicated by the presence of data points 
up to, and even beyond 0.1 THz, as well as the wider 
probability distribution observed at low frequencies. In 
order to study the convergence of the lifetimes towards 
a single value, we group the data for different A and the 
same a together, and calculate the standard deviation 
cr{5ijj) of the resonances width Sijj. Figure |4f illustrates 
the results of this analysis. In perfect agreement with 
our theoretical predictions, we observed a significant nar- 
rowing of the linewidth distribution above the threshold 
for chaos. (Fig. |4ti). It is worthwhile highlighting that 
this convergence of the linewidth towards a single value 
does not depend on A, but only on a, which is a clear 
experimental demonstration that the phenomenon relies 
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FIG. 4. Summary of 2D experimental results: (a)-(c) SEM micrographs of the sample geometry for (a) a = 0, (b) a = 0.75 
and (c) a = 1; (d) experimental spectrum (circle markers) and theoretical reconstruction via wavelet multiscale analysis (solid 
line) for A = 400/im^ and a — I; (e) probability distribution P{Suj) of the resonances width Su calculated for the fully chaotic 
resonator a = 1; (f) standard deviation a{Suj) of the resonance widths versus a. 
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FIG. 5. Summary of the 3D experimental results with deformed microspheres: (a) sketch of mono-layer sample fabrication by 
convective self-assembly; (b)-(c) furnace for sample heating and deformation with mechanical pressure; (d)-(e) SEM micrograph 
of a microsphere in the original (d) and deformed (e) case; (f) average normalized absorption of the deformed microsphere 
< Ad >=< Ad/A >, measured for different wavelengths and normalized with respect to the undeformed case A. 



entirely on the chaotic properties of the motion of light. 



3D experiments with deformed microspheres: 
chaos-enhanced broadband absorption of light 

A true comparison between model and experiment re- 
quires a direct measurement of the energy in the system. 
We accomplished this final step in a three dimensional 
geometry, which incidentally also proves how the physics 



of chaotic resonators is independent of the dimension- 
ality of the problem. We employed polystyrene micro- 
spheres, and studied their absorption when their spher- 
ical shape gets deformed. Our sample consists of a low 
density monolayer of microspheres (Fig. [5^), which are 
being deformed by mechanical compression (Fig. [Sj)) in 
order to obtain an asymmetric shape such as that shown 
in Fig. [T}d. The monolayer fabrication was performed by 
employing convective assembly [23, 24], where the micro- 
spheres are self- assembled on a substrate via a deposition 
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blade (Fig. [5^). The blade height was set at 12/im from 
the substrate (the sphere size is « lO/xm) to ensure the 
formation of a single layer. Following deposition, a glass 
slide was placed on top of the microspheres for heating 
and applying mechanical pressure. The deformation was 
realized by heating the system slightly beyond the glass 
transition temperature Tg of polystyrene {Tg = 94° Cel- 
sius [Isj), which softens the microspheres sufficiently to 
deform their shape (Fig. [Sjj-c). Compression and heat- 
ing were performed with an Instron 5960 dual column 
tabletop universal testing system. The pressure force was 
ramped to SOON over 2.5 mins. Figures [5ji-e show opti- 
cal micrographs of a microsphere before (d) and after (e) 
compression, highlighting the asymmetric deformation. 
We evaluated the energy harvesting capacity of this sys- 
tem by performing absorption measurements on a single 
sphere in both deformed and undeformed conditions. In 
order to acquire a sufficiently large statistics from dif- 
ferent input conditions, we placed our sample on a go- 
niometric stage and, through a series of pump and probe 
measurements, we calculated the absorption of the sphere 
at various illumination angles in the range of ±30 de- 
grees. During each measurement, we illuminated the 
sphere with a broadband source (with bandwidth « l/xm 
centered at 700nm) and measured the absorption A from 
the relation A — 1 — T — being T, R transmittance 
and reflectance of the microsphere, respectively. To prop- 
erly collect the scattered light, transmission and reflec- 
tion spectra were measured in the near field with a lens 
with a high numerical aperture. The data was acquired 
with an Ocean Optics QE65000-FL spectrometer. Fig- 
ure [5f shows the average absorption < Ad >=< Ad/ A > 
of the deformed system, being Ad and A the absorption 
measured in the deformed and undeformed case, respec- 
tively. For every angle we obtained an increased absorp- 
tion due to deformation, with a variance within a few 
percent. This originates from the larger electromagnetic 
energy stored by the microsphere in the deformed case, 
which led to a higher absorption in the entire absorption 
window of polystyrene (Fig. [5f). Near the wavelength 
A = 450 nm, in particular, the average absorption in- 
creases by « 12%, while in the region where the optical 
spectrum contains the maximum power (« GOOtito), the 
absorption shows a broadband increase of a approxima- 
tively 6 — 8%. Outwith the absorption frequency win- 
dow of polystyrene (Fig. [5] solid lines), as expected, the 
chaotic energy accumulation does not induce any mea- 
surable variation of the absorption. 



Discussion 

Our analytical, numerical and experimental results ad- 
dressed the problem of light trapping in chaotically de- 
formed resonators and demonstrate the possibility of ex- 
ploiting chaos as a new avenue for energy harvesting. 
Our work began with a single numerical ab-initio ex- 
periment, where we demonstrated broadband energy ac- 



cumulation in an optical resonator once its shape was 
deformed to support the chaotic motion of light. We 
then developed an analytical approach based on a rigor- 
ous formulation of time dependent coupled mode theory, 
which highlights the principal phenomena of chaotic en- 
ergy trapping. This approach allowed us to unveil the 
presence of the fundamental thermodynamic principle of 
energy equipartition. Equipartition, in our case, is man- 
ifested by the uniform distribution of energy among all 
degrees of freedom of the chaotic resonator, i.e. the cav- 
ity modes, and is evident from the convergence of the 
lifetime of the modes towards a single value. A single life- 
time means that each mode can store the same amount of 
energy, which overall leads to a significant energy buildup 
in the resonator. The process is equivalent to the Brow- 
nian motion of particles in a liquid, where each particle 
carries the same amount of energy. The analogy goes 
even further. In a liquid, the particles always reach a 
uniform distribution regardless of the shape of the ves- 
sel they are contained in. In our optical analogue, the 
steady state of chaotic modes also does not depend on 
the particular realization of the resonator, but only on 
its 'macroscopic' geometry, whose only requirement is to 
support chaotic trajectories for the trapped light. A large 
series of parallel ab-initio simulations perfectly confirmed 
our theoretical predictions. 

We then conducted an experimental campaign to demon- 
strate both the main physical principles and the possible 
applications of this concept. In particular, we started 
by employing two-dimensional (2D) photonic crystals 
and realized stadium-shaped resonators with different 
degrees of chaoticity. Using pump and probe experi- 
ments, we demonstrated the convergence of the modal 
decay times towards a single value, which is at the ba- 
sis of energy equipartition and chaotic energy harvest- 
ing. We then considered the realization of a practical 
three-dimensional (3D) geometry where we demonstrated 
enhanced light trapping through shape deformation di- 
rectly. We employed polystyrene microspheres and al- 
tered their shape through heating and mechanical com- 
pression. We then measured the absorption of the sys- 
tem, and experimentally demonstrated a broadband in- 
crease across the entire absorption window of the mate- 
rial by deforming the spheres form spherical to chaotic. 
The measured absorption increase was as high as 12%, 
which is quite remarkable considering the low absorption 
of polystyrene and the low refractive index of the mate- 
rial which only affords weak light trapping . 
Besides the obvious implications at the fundamental 
level, where we demonstrate the existence of a funda- 
mental principle of thermodynamics in the framework of 
Photonics, our results also have real- world practical im- 
plications. The cost of many semiconductor devices, e.g. 
LEDs and solar cells, is determined to a significant extent 
by the cost of the material. We show that the function- 
ality of a given geometry, here exemplified by the energy 
that can be trapped in the system, can be enhanced up 
to six-fold by changing the shape alone, i.e. without in- 
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FIG. 6. (a) light-cavity interaction TDCMT scheme: S{t) is 
the input source, R{t), T{t) are reflection and transmission 
signals, respectively, at is the fc-th mode in the resonator 
and 7 is the coupling coefficient between the resonator and 
the external environment; (b) normalized entropy K versus 
chaotic degree a for the resonator defined by Eq. ((!]). 



creasing the amount of material and without increasing 
the material costs. Furthermore, a chaotic system is eas- 
ier to fabricate as the tolerances are relaxed; the only 
requirement is that the resonator supports chaotic light- 
paths, which is much easier to achieve in practise than 
manufacturing a resonator which is perfectly regular. 



modes obey the following evolution equations [18|: 



dak 
~dt 







iujk - 









ak+i-S{t), fce[l,...,n], (4) 



with "f = ^ [l^, while reflection R = J2k^k and 

transmission T — T^. are given by the following ex- 
pressions: 



Rk 



1 

— ak, 



Tk 



-+S, 



(5) 



with Rk and Tk the reflection and transmission of the 
fc— th mode, respectively [11]. In the presence of a single 
frequency excitation S = e*'^*0(i) switched on at i = 0, 
with 9(t) being the Heaviside function, Eqs. Q are 
readily solved for each ak and read: 



ak 



1 e*' 



e i{uj - Wfe) + — 



(6) 



For a broadband source, S — J (iwe'"*8(t), the total elec- 
tromagnetic energy £ = jafep stored in the cavity is 
readily found to be: 



Methods 

Ab-initio computations. — Numerical simulations 
were performed with our NANOCPP Maxwell solver 
(www.primalight.org). NANOCPP is a two/three- 
dimensional FDTD Maxwell code based on a paral- 
lel implementation of the Yee FDTD algorithm. The 
code employs uniaxial perfectly matched layers (UPML) 
to efficiently absorb outgoing waves and provides both 
hard and planewave sources [the latter via total field 
scattered field (TFSF) formulation]. NANOCPP effi- 
ciently scales over hundreds of thousands of processors, 
and is specifically designed for massively parallel elec- 
tromagnetic computations. Our simulations consisted of 
a total of 20 Million of single-cpu computational hours, 
with each simulation employing 4096-8192 computational 
cores on the KAUST Shaheen computational cluster and 
a resolution of at least 20 points per internal wavelength. 

Time dependent coupled mode equations. — Figure 
shows a TDCMT model of the light-resonator sys- 
tem. In this scheme, the environment is represented as a 
waveguide side coupled to the resonator. The input sig- 
nal S{t) interacts with the resonator through the coupling 
coefficient 7 and transfers energy to the modes of the cav- 
ity ak{t) = Afce'"^"*"*/^" (fc G [l,iV]), being ujk the fre- 
quency, the lifetime, Ak the amplitude and £k = lofcp 
the energy carried by the fc— th mode. The total energy 
stored in the resonator is £ = X^fc^fe- Decay rates 
can be further decomposed as follows — = — + — , with 

— the intrinsic decay constant of the mode and — the 
decay constant due to coupling with the source. Cavity 



H = 



duj 1 

'''^ k 



2 cos[{uJk ~ uj)t]t 



(7) 



Equation ([7]) can be further simplified as the integral 
yields significant contributions only for uj ujk, and we 
obtain: 



(8) 



which represents the first of Eq. ([2]). The power Vk 
transferred into the fc— th mode is conversely evaluated 
from the energy balance equation: 



d\ak\ 
dt 



Vk - \S[' 



\Rk? 



\Tk\\ 



(9) 



that, in conjunction with Eqs. (IS])-®, yields the second 
of Eqs. ©. 

Characterization of chaos. — We quantitatively de- 
scribe the chaoticity of light motion by evaluating the rel- 
ative area of the resonators phase space that encompasses 
chaos. This is achieved by first calculating the distribu- 
tion of the Lyapunov exponent [2^ in phase space, and 
then performing a weighted summation by assigning T' 
if the Lyapunov exponent is positive, and '0' otherwise. 
The resulting quantity K can be regarded as a normal- 
ized version of the Kolmogorov-Sinai entropy [26l.[27|. In 
particular, when K — Q, the resonator dynamics exhibits 
no chaos and the resulting motion is totally reversible, 
while for K = 1 the phase space is totally chaotic and 
all input conditions lead to chaos. Values of K between 
these two limiting conditions indicate a phase space par- 
tially chaotic, with K measuring the relative area of the 
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chaotic sea with respect to the reversible portion of the 
dynamics. Figure ^jp displays the behavior of K for the 
billiard of Eq. [T] For a lower than the threshold value 
a* — 0.75, no chaos is observed in the structure, while for 
a > a* strong chaos is generated through the shape de- 
formation. At a = 1.5, the structure is fully chaotic and 
the phase-space is totally dominated by a single chaotic 
sea. 
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